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Abstract. We report on the structure of the nuclear star cluster in the 
innermost 0.16 pc of the Galaxy as measured by the number density profile of 
late-type giants. Using laser guide star adaptive optics in conjunction with the 
integral field spectrograph, OSIRIS, at the Keck II telescope, we are able to 
differentiate between the older, late-type (~ 1 Gyr) stars, which are presumed 
to be dynamically relaxed, and the unrelaxed young (~ 6 Myr) population. This 
distinction is crucial for testing models of stellar cusp formation in the vicinity 
of a black hole, as the models assume that the cusp stars are in dynamical 
equilibrium in the black hole potential. In the survey region, we classified 77 
stars as early-type and 79 stars as late-type. We find that contamination from 
young stars is significant, with more than twice as many young stars as old stars 
in our sensitivity range (K < 15.5) within the central arcsecond. Based on the 
late-type stars alone, the surface stellar number density profile, oc R~ r , is 

fiat, with T = — 0.26±0.24. Monte Carlo simulations of the possible de-projected 
volume density profile, n(r) cx r~ 7 , show that 7 is less than 1.0 at the 99.7 % 
confidence level. These results are consistent with the nuclear star cluster having 
no cusp, with a core profile that is significantly flatter than predicted by most 
cusp formation theories, and even allows for the presence of a central hole in the 
stellar distribution. Here, we also review the methods for further constraining 
the true three-dimensional radial profile using kinematic measurements. Precise 
acceleration measurements in the plane of the sky as well as along the line of 
sight has the potential to directly measure the density profile to establish whether 
there is a "hole" in the distribution of late- type stars in the inner 0.1 pc. 



1. Introduction 

Over 30 years ago, theoretical work suggested that the steady state distribution 
of stars may be significantly differ ent for clusters with a massive black hole at 
the center th an those without (e.g. iBahcall fc WolJll976l ; ICohn fc Kulsrudll 19781 ; 



Youngl 1 19801 ) . Stars with orbits that bring them within the tidal radius of the 



black hole are destroyed and their energy is transferred to the stellar cluster. In 
the steady state, this energy input must be balanced by the contraction of the 
cluster core, which appears as a steeply rising radial profile in the number density 
of stars toward the cluster center. This radial profile is usually characterized by a 
power law of the form n(r) oc r~ 7 , with a power law slope, 7, that is steeper than 
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that o f a flat isothermal core. For a single-mass stellar cluster, iBahcall &: Woli 
( 19761 ) determined the dynamically relaxed cusp will have 7 = 7/4. The presence 
of such a steep core profile, or cusp, is important observationally because it may 
represent a simple test for black holes in stellar systems where dynamical mass 
estimates are difficult, such as in the cores of galaxies. The stellar cusp is 
also important theoretically as it is a probable source of fuel for the growth 
of supermassive black holes and its presence is of ten assumed i n simulations of 
stellar clusters having a central black hole (e.g., iMerritt 20051 . and references 

therein). 

Theoretical work has progressed from the early simulations bv lBahcall &: W olf 

(| 19761 ) to include many complicating effects, including multiple masses, mass 
segregation, and stellar collisions, on the d ensity profile of stellar clusters in the 
presence of a supermassive black h ole (e.g.. IBahcall k. WolJll977l ; iMurphv et al 



1991 ; Alexander fe Hopmanl 12009). These theories predict cusp slopes within a 



black hole's gravitational sphere of influence ranging from 7/4 > 7 > 7/3 for a 
multiple mass population to as shallow as 7 = 1/2 for a collisionally dominated 
cluster core. An important assumption of all cusp formation models is that the 
stellar cluster be dynamically relaxed. Without this assumption, the stellar dis- 
tribution would show traces of the cluster origin in addition to the influence of 
the black hole. 

The Galactic center is an ideal place to test these theories of cusp formation 
as it contains the nearest exam ple of a superma s sive black hole (Sgr A *, with a 
massM. = 4.1±0.6x 10 6 M ? . lGhez et~alll2008l : [Gillessen et al.ll2009l ). Located 
at a distance of only 8 kpc, the radius of the sphere of influence of the Galactic 
center black hole (~ 1 pc) has an angular scale of ~ 25" in the plane of the sky, 
two orders of magnitude larger than any other supermassive black hole. In order 
to plausibly test stellar cusp formation theories, the stellar population used to 
trace the cusp profile must be older than the relax ation time, which is of order 1 
Gyr within the sphere of influence of Sgr A* (e.g., iHopman Sz Alexander| [2006). 
In the Galactic center, late- type red giants (K to M) are the most promising trac- 
ers of the stellar distribution since they are > 1 Gyr old, bright in the NIR, and 
abundant. On parsec size scales, these stars dominate the flux, so early seeing 
limited observations of the surface brightness profile in the NIR at the Galac- 
tic center were successfully used to show that the struc ture of the nuclear star 
clust er at large scales has a density power law of 7 « 2 (jBecklin Neugebauer 
1968). Integrated light spectroscopy of the CO absorption band-head at 2.4 /j,m 
(which is dominated by red gi ants) also confirmed this slope down t o within 
0.5 pc from the black hole (|McGinn et all 1 19891 : lHaller et all I1996T). These 



measurements, as well as integrated light spectroscopy from Figer et al. (2000) 



of the inner 0".3, found a lack of CO absorption in the inner ~ 0.5 pc of the 
Galaxy. It was unclear at the time whether this represented a change in the 
stellar population with fewer late-type stars in the inner region or a change 
in the stellar density profile. Because the integrated light spectrum is biased 
toward the brightest stars, contamination from a few young Wolf-Rayet (WR) 
stars (age < 6 Myr), such as the IRS 16 sour ces identified in the early 1990s 
locat ed between 1-2" from the black hole (e.g., Allen et al.lll990l : iKrabbe et al.l 



can significantly impact measurements of the underlying stellar density. 
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Here, we report on an update of our high angular resolution spectroscopic 
survey of the central 0.15 parsec of the Galaxy using laser guide star adaptive 
optics with the OSIRIS IFU at the Keck II telescope. The survey reaches a 
completeness of 40% at K' ~ 15.5, two magnitudes fainter than previously 
reported spectroscopic surveys of late-type stars in this region. We find that the 
surface number density of late-type stars is sign ificantly flatter than the range of 
power laws predicted bv lBahcall fc Woll ( 19771 ) (Section [3]]) • This measurement 



rules out all values of 7 > 1.0 at a confidence level of 99.7%. In Section d3 we 
discuss possible dynamical effects that may lead to the flat observed slope and 
implications for future measurements of the stellar cusp at the Galactic center. 

2. Observations 

Near-IR spectra of the central 4" of the Galaxy were o btained between 2 006 
and 2009 using the OSIRIS integral field spectrograph ( Larkin et al.l 120061 ) in 



conjunction with t he laser guide star adaptive optics (LGS AO ) system on the 
Keck II telescope ( Wizinowich et al. 20061 ; Ivan Dam et al.ll2006l ). Details of the 



observations are given in (IDo et al.ll2009r). In this paper, we include two addi- 
tional fields not reported in ()Do et al.ll2009l ). the northwest and southwest field, 
extending complete azimuthal coverage in the survey out to ~ 3" from Sgr A* 
(see Figure [I]). The two additional 2" x 3" fields are centered at distance of 
-1.99", 2.4" (northwest) and -2.99", -1.12" (southwest) from Sgr A* (obtained 
in 2009 May). 

3. Results 

With spectroscopic identification of the stars brighter than K' ~ 15.5, we are 
able to separate the presumably dynamically relaxed old stars from the unre- 
laxed young population. The stellar number density profile is constructed in 
radial bins of 0."5 over the area that was sampled in this survey out to a radius 
of ~ 3", with error bars scaling as y/~N, where N is the number of stars in each 
bin. It is clear that, while the number density of early-type stars increases quite 
rapidly toward the central arcsecond, the late-type stars have a very flat distri- 
bution (Figure [Q). Outside the central arcsecond, the projected number density 
of early-type stars drops to about half that of the late-type stars. 

Potential systematic errors to the radial profile measurement may arise 
from variable extinction and incompleteness between each field. In order to 
study how these two factors impact our result, the radial profile was computed 
using stars with extinction corrected magnitudes only over regio ns that are at 



least 30% complete. To do so, the extinction map reported by (jSchodel et al.l 
l2009h was used to correct for the variations in extinction between each star 
by adding Ak — 3.0 to each star to bring them to the same canonical Ax = 
3.0 extinction to the Galactic center. Then, the K' luminosity function was 
recomputed for each field to determine the K' magnitude bin beyond which the 
completeness falls to less than 30%. This separates the fields into three groups 
in terms of their 30% completeness : all fields are at least 30% complete down 
to K' = 14.5; all fields except W, NW, and SW are complete to K' = 15.5; 
the central field is complete to K' = 16.0. We calculate the radial profiles in 
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Figure 1. The currently surveyed region overlaid on a K' (2.2 /mi) image 
of the Galactic center taken in 2007. Sgr A* is marked at the center with 
a *. Spectroscopically identified early (blue triangles) and late- type (red cir- 
cles) stars are marked. Each field is enclosed by dotted lines. Some spectral 
identifications are outside of the marked lines because they were found at the 
edge of the dithers. 

each magnitude bin using only fields that are at least 30% complete at that 
magnitude. For all stars with K' < 14.5, we use all the fields, while for stars 
with 14.5 < K' < 15.5, the western field was dropped from the measurement; 
the 15.5 < K' < 16.0 magnitude bin was not included because only the central 
field is complete to at least 30% at that magnitude. The radial profiles from the 
K' < 14.5 and 14.5 < K' < 15.5 magnitude bins were then summed to produce 
the final measurement of the surface density, with errors added in quadrature 
(see Figure [2]). The early-type stars have T young = —0.90 ± 0.27, while for the 
late-type stars, r o ^ = —0.26 ± 0.24. Neither of these values differ significantly 
from the case without extinction and completeness correction. The rest of the 
analysis will use these values for the measured surface number density profile. 

4. Discussion 

The effects of projection of a three-dimensional cluster onto the plane of the sky 
must be taken into account for quantitative interpretations of the measured sur- 
face density profile. The projection onto the sky plane of a spherically symmetric 
cluster can be done using the integral: 




(1) 
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Figure 2. Plot of the surface number density as a function of projected 
distance from Sgr A* in the plane of the sky for different populations: old 
(late- type, red), young (early-type, blue), and total number counts from K' 
imaging. These radial profiles have be en corrected fo r completeness and ex- 
tinction using the method detailed in (|Do et al.ll2009|) . 



where R is the distance from the center of the cluster in the plane of the sky, r is 
the physical radius of the star cluster, E(-R) is the projected surface density, and 
n(r) is the volume number density. Because this integral diverges for all 7 < 1.0, 
where n(r) oc r -7 , the density in the outer region of the cluster must fall off 
steeply to accommodate a shallow power law at the center of the cluster. While 
we do not observe a break in the cluster density profile within our 4" survey, 
one was observed by ISchodel et all ( 20071 ) and others at scales larger than our 



survey. In order to account for the outer region of the star cluster, we therefore 
model the number density profile as a broken power law: 



Using the measurements from Schodel et al. ( 20071 ). we set ri, rea k = 8" and 72 



1.8 (the deprojected values of Rbreak = 6" and T = 0.8, from the surface density 
measurements using this broken power law model) to constrain the parameters 
of the outer cluster density profile. We assumed that contamination from young 
stars is less severe in the region outside of our survey, so using an outer cluster 
density profile measured from total number counts should introduce only a small 
bias. For example, the density of bright yo ung stars fall s off relatively steeply as 
R- 2 (iPaumard et al.ll2006l : iLu et alJ Eo09. but see also iBuchholz et al.1 (|2009l )). 



To determine the constraints on the physical density profile in the inner region, 
we performed Monte Carlo simulations of clusters with 71 between -2.0 and 2.0. 
For each value of 71, 10 4 realizations were performed. In each realization, 74 
stars (the number of late-type stars used to measure the radial profile) were 
drawn from the broken power law distribution after correcting for the limited 
area observed in this survey. The locations of these stars were then binned and 
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fitted the same way as the data. We find that for —2.0 < 71 < 0.5, the projected 
inner radial profiles are flat with T ~ 0. Using the relationship between 71 and 
r as determined from the MC simulations and the observed inner radial profile, 
we can set an upper limit of 1.0 to the value of 71 at the 99.7% confidence level. 



4.1. Mechanisms for cusp depletion 

Although the slope measurement in this survey cannot constrain whether there 
is a 'hole' in the distribution of late-type stars or just a very shallow power law 
within the central 4", the inferred slope is significantly flatter than the range 
of 71 between 7/4 and 3/2 predicted by iBahcall & Wolfl (|l977l ). In fact, a flat 
core density profile without a cusp can fit the dat a equally well. This r esult 
is consistent with si milar recent measurements from iBuchholz et al.l ( 20091 ) and 
Bartko et al.l ( 20101 ). Below, we discuss several plausible dynamical scenarios 
that may deplete the number density of the late-type giants observed in this 
study, as well as possible observational prospects for making further progress. 



Mass segregation The true population of dark mass that consists of stellar 
remnants such as neutron stars or stellar black holes is unknown, but theoreti- 
cally, there i s a strong cas e for a large population of such remnants due to mass 
segregation ( Morris! 1993). As these dark remnants mig rate inward, they will 
scatte r the lighter stars outward. Recent simulations bv lAlexander fe Hopmaii 
(2009) showed that a population of ~ 10 M stellar remnants can sink into 
a much steeper density distribution (2 < 7 < 11/4) than the lighter stars 
(3/2 < 7 < 7/2), if the stellar remnants are relatively rare compared to stars. 
This results in a very dense cluster of stellar remnants, but the prediction of 
the stellar cusp from stron g mass segregation does not differ substantially from 
the IBahcall fc Wol (|l977l 1 values. Thus, mass segregation cannot be the only 
mechanism responsible for the inferred lack of a cusp in the evolved red giant 
population. 



Envelope destruction by stellar collisions The M and K type giants in this 
survey have radii that vary between 20 to 200 R©, which results in substantial 
cross sections for collis i ons with other st ars in this high density environment (e.g. 
Bailev fc Davieslll999l : foale et al.1 l2009h . This may lead to a preferential deple- 
tion of red giants toward the center, which would result in a biased measurement 
of the underlying stellar cusp slope, since most of the surviving stars in the cusp 
would be unobserved main sequence stars with smaller radii. The depletion of 
giants has been suggested in the past based upon low spatial re solution spec- 
troscopy of the CO b and-head at 2.29 // m dMcGinn et al.lll 989) as well as m 
narrow band imaging (|Genzel et al.ll2003l : IBuchholz et al.l 12009). Theoretically, 
a collisionally dominated cusp can be as shallow as 7 = 0.5, which is within the 
range of slopes allowed by our upper limits, for an isotropic distribution, (i.e., 
for a distribution function that depends only on energy), although this result no 



longe r holds in the case of a more general distribution function (jMurphv et al 
1991). While the collisional destruction of giants can operate efficiently within 
the region observed in this survey, the effectiveness of this mechanism drop s 
very quickly with distance from the black hole. For example, iDale et ail ( 2009) 
found that collisional destruction of giants becomes less likely beyond about 0.1 
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pc from Sgr A*. While this survey covers a region ou t to ~ 0.16 pc i n pro - 
jection from Sgr A*, previous observations such as from iBuchholz et al. ( 20091 ) 
suggest that the flattening in the late-type giant radial profile extends out even 
further to ~ 0.24 pc in projection. If the assumptions about the steep distribu- 
tio ns of stella r mass black holes and other stellar remnants in the simulations 
of iDale et all (|2009l ) are correct, then this would suggest that the collisional 
destruction of giants is not the dominant mechanism for clearing the cusp of 
stars. Observational constraints on the effectiveness of this mechanism can be 
made by comparing the radial density profiles of giants of varying stellar radii. 
Collisional destruction should preferentially remove stars having larger radii at 
a given distance from the black hole. 



IMBH or binary black hole A more dra matic change to the stell ar cusp can be 
achieved through the infall of an IMBH. Baumgardt et al. ( 20061 ) performed N 
body simulations of the effects on the central star cluster with an initial cusp 
slope of 7 = 7/4 due to the infall of IMBHs of mass 3 x 10 3 M and 10 4 M . 
They find that the IMBH will deplete the cusp of stars, which will cause the inner 
portion of the cluster to resemble a core profile after the IMBH has spiraled into 
the central black hole, and that it would take over ~ 100 Myr after the IMBH 
infall to replenish the stellar cusp. For the case of a 10 4 M IMBH, the density 
profile has a flat core, with a core radius of ~ 0.2 pc, which is r oughly consistent 
with that observed in this survey and by IBuchholz et~all (120091 ) . 

The existence of a binary massive black hole in the Galactic center can 
also contribute to the loss of stars in this region. A star passing within the 
orbit of the companion massive black hole will undergo dynamical interac- 
tion s with the binary, which can result in the star achieving ejection veloc- 
ity ( Merritt &: Milosavlievic 20051 ). In addition, similar to that of an infalling 
IMBH, a merger with a massive black hole system more recently than the dy- 
namical r elaxation time will r esult in a system that is still out of equilibrium at 
this time (jMerritt et al.ll2007l ). 

Current measurements impose only modest limits on th e possibility of an 
infalling IMBH and a binary black hole. iTrippe et al. ( 20081 ) found no obvious 
kinematic signatures of a recent infall of an IMBH using 3D velocity measure- 
ments of the late-type stars, however, a more complete study of the phase space 
distribution of the stars will b e necessary to ascertain the l ikelihood of an IMBH 
clearing the cusp of stars (see iGualandris Sz Merrittl ( 20091 ) for a more extensive 
discussion of current constraints on the existence of an IMBH at the Galactic 
center). Current measurements of the reflex motion of the black hole at radio 
wavelengths place the limit to the mass of a c ompanion black hole at < 10 4 M 
with semi-major axes 10 3 AU < a < 10 5 AU (|Reid fc Brunthaledl2004T ). Precise 
measurements of the possible reflex motion of the black hole using stellar orbits 
can provid e stronger constra ints in the future on the presence of a companion 
black hole (|Ghez et al.ll2008l ). 



4.2. Constraining line of sight distance with kinematic measure- 
ments 



The degeneracy in the projection from a three-dimensional cluster density profile 
onto the plane of the sky must be overcome in order to better constrain the 
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theories for the depletion of the stellar cusp. This degeneracy can be resolved 
by measurements of the line of sight distance to each star. The star cluster in 
the Galactic center has the unique advantage of having very precise kinematic 
measurements, which with some assumptions about the potential, can be used 
to derive the line of sight distances of the stars, z, from the black hole. 

The simplest method of determining the limits to z for each star comes 
from comparing the total velocity, Vtot of the star to the escape velocity at 
the star's projected distance, R from the black hole. If we assume the major 
contributor of the potential to be the black hole mass, Mbh, and that the stars 
are gravitationally bound, then: 



For the stars at the center of the nuclear star cluster, the line of sight 
distance can be measured directly using accelerations. With the assumption 
that the black hole mass dominates the potential, the true radial distance, r, 
from the black hole can be simply computed using Newton's law of gravity: 

GM BH 

a = —pi — v 4 ) 

where a is the magnitude of the acceleration of the star. This assumes that the 
acceleration can be measured in both the plane of the sky and along the line 
of sight. It is useful to decompose these two components, because accelerations 
in the plane of the sky are measured using astrometry, while the accelerations 
along the line of sight are measured spectroscopically with different measurement 
uncertainties. The line of sight distance can be computed from the plane of the 
sky acceleration, 02/5, using: 

„, n _ GM BH p 

a2D - -(#+32)3/2 « , , 

Z = (( GM BH R }2/3 ^2)1/2 W 

With the measurements of a,2D alone, ±z is degenerate (the star can be in front 
or behind the black hole), but is not of concern for measuring the radial profile 
of the number density. More problematic however is that as R2D approaches 0, 
OL2D also approaches (as the acceleration vector always points toward the black 
hole) , resulting in less sensitivity to the measurement of z for a given astrometric 
precision. In contrast, acceleration along the line of sight from spectroscopic 
measurements are more sensitive for stars close to the black hole in the plane of 
the sky: 

a z = ~ {R 2 +z B 2ys/2 z (6) 

However, z is degenerate for measurements of a z over a range of z. For example, 
if a z is small, z can be very close to zero, or z can be very large. Limits on a,2D are 
required to break this degeneracy. Figure [3] summarizes our sensitivity to z given 
our average acceleration uncertainties in a z and a,2D- Currently, we should be 
able to detect all stars within 1" in volume around Sgr A*. Beyond this volume, 
we have some sensitivity out to ~ 2.5". With a longer time baseline and lower 



9 



3.5 



3.0 



2.5 



2.0 
o 

1 1.5 

N 



1.0 
0.5 



0.0 



a2d sensitivity 
az sensitivity 



I i I . : I 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
R2D distance (") 



Figure 3. The sensitivity to the line of sight distance for stars based on 
the average acceleration errors along the line of sight (a 2 D) and in the plane 
of the sky (a z ) are shown in the shaded regions. The dotted lines show the 
volume of space encompassed within 1", 2", and 3" from Sgr A*. 



acceleration uncertainties in measurements of a,2D, we are able to measure the 
z distances of stars better than with a z in the survey region. Further details on 
the measurements and the three dimensional profile of the cluster will be given 
in Do et al. (in prep). 



5. Conclusions 

We completed a survey using high angular resolution integral field spectroscopy 
of the inner ~ 0.16 pc of the Galaxy in order to measure the radial profile of 
the late-type stars in this region. The survey reached a depth of K' < 15.5 and 
is able to differentiate between early and late-type stars. We find a late-type 
stellar surface density power-law exponent V = —0.26 ± 0.24, which limits the 
volume number density profile slope 7 to be less than 1.0 at the 99.7% confidence 
level, and even allows for the presence of a central drop in the density of late- 
type giants. Given the current measurements, we cannot yet determine whether 
the distribution of observed giants is representative of the distribution of stellar 
mass. Being able to infer the underlying dynamically relaxed stellar population 
will be crucial in order to establish whether the Galactic center is lacking the type 
of stellar cusp long predicted by theory. Obtaining an unbiased measurement of 
the stellar distribution is important because such cusps have considerable impact 
on the growth of massive black holes as well as on the evolution of nuclear 
star clusters. Progress in achieving this goal will be possible with improved 
kinematics and spectral coverage in order to break the degeneracy in the surface 
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number density profile to better establish the three-dimensional distribution of 
the stellar cluster. 
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